Carbon clusters: from ring structures to nano-graphene 
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The lowest energy configurations of Cn{n < 55) clusters are obtained using the energy minimiza- 
tion technique with the conjugate gradient (CG) method where a modified Brenner potential is 
invoked for the carbon and hydrocarbon interaction. We found that the ground state configuration 
consists of a single ring for small number of C atoms and multi-ring structures are found with 
increasing n, which can be in planar, bowl-like or cap-like form. Contrary to previous predictions, 
the binding energy Eh does not show even-odd oscillations and only small jumps are founded in 
the Eh{n) curve as a consequence of specific type of edges or equivalently the number of secondary 
atoms. We found that hydrogenation of the edge atoms may change the ground state configuration 
of the nanocluster. In both cases we determined the magic clusters. Special attention is paid to 
trigonal and hexagonal shaped carbon clusters and to clusters having a graphene-like configuration. 
Trigonal clusters are never the ground state, while hexagonal shaped clusters are only the ground 
state when they have zigzag edges. 

PACS numbers: 71.10.Pm, 73.21.-b, 81.05.Uw 



I. INTRODUCTION 

The fabrication of single layers of graphite - known as 
graphene^ has resulted in an intense interest in the struc- 
ture and properties of stable graphene-like structures. 
Questions concerning the formation mechanisms and the 
possible stable structures, especially the most stable 
ones, has focused attention on small carbon clusters both 
by experimentalists and theorists. Even with modern 
probe techniques^ it is often difficult to obtain unambigu- 
ous structural information from experiment^ ^ for small 
size clusters. On the theoretical side, the study of carbon 
clusters can be divided into two main categories: calcu- 
lations using abinitio techniquesSHSHni a nd calcula tions 
based on empirical interatomic potentials^^'^^l^^^. Ac- 
curate ab-initio calculations performed by Raghavachari 
and Binkley (Ref. 7) obtained the structure and energy 
of small carbon clusters Cn{n = 2 — 10) and predicted 
an odd-even alternation in the nature of the cluster ge- 
ometries with the odd-numbered clusters having a linear 
structure and the even-numbered clusters preferring an 
irregular cyclic structure. 

There are two kinds of empirical interatomic carbon 
potential functions that are widely used. One is the em- 
pirical potential proposed by Brenner, which was based 
on the Abell-Tersoff bond-order expression and fitted on 
diamond, graphiti crystals and small hydro carbonsl. 
The validity and reliability of the Brenner potential 
to the study of fullerene, nanotubes and carbon clus- 
ter formati on h ave been demonstrated in a range of 
applications^^EMU] With such potential function C. 
Zhang et al. (Ref. 13) used a genetic algorithm (GA) to- 
gether with a simulated annealing (SA) method to study 
the geometry ofCn{n = 2 — 30) clusters. D. H. Robertson 
et al.^ used molecular dynamics (MD) to obtain three- 
ring structures up to fullerenes by modifying one of the 
parameters in the Brenner potential function. Y. Yam- 
aguchi and S. MaruyamaP^, and S. Maruyama and Y. 



Yamaguchi^^ used molecular dynamics (MD) to study 
the formation of fullerenes out of 500 carbon atoms that 
were generated randomly in a cubic box by omitting one 
of the parameters in such a Brenner- type potential func- 
tion. Y. Wang et aL^^used a time-going backward quasi- 
dynamics method with the Brenner potential in order to 
optimize cluster structures. K. Mylvaganam et al.^^ ap- 
plied MD to analyze carbon nanotubes and their mechan- 
ical properties. W. Cai et al.^^ used a global optimization 
algorithm to study structural properties of Cn(n = 2—71) 
clusters. 

Another kind of potential function is the semi- 
empirical tight- binding (TB) potential which was based 
on the extended Hiickel molecular orbital approximation. 
Such kind of potentials were used to obtain equilibrium 
geometries for carbon clusters of arbitrary size^^ ^^. 

Recently, magnetic order in different shaped graphene 
clusters, that arise from size effects or from topological 
frustration, has been proposed as building elements for 
logic gates in novel ultrafast high-density spintronic de- 
vices that can work at room temperature^^. In these 
calculations the shape of the clusters and ordening of 
the carbon atoms in the clusters were assumed to be the 
same as in the ideal graphene sheet without relaxation 
of the C-atom positions and without questioning if such 
clusters can be stabilized or not. 

Here, using energy minimization with the conjugate 
gradient (CG) method, we obtain the stable states of 
different size clusters and obtain the lowest energy config- 
uration of the carbon atoms for quasi-2D arrangements. 
We investigated the equilibrium configuration of different 
size clusters and examined when graphene lattice-type 
structures can be formed. Both even- and odd-numbered 
neutral carbon clusters Cn were systematically studied in 
our previous work^^ for(n = 2 — 10). All possible stable 
configurations were found, and several new isomers were 
predicted. It should be pointed out that the dependence 
of the binding energy for linear and cyclic clusters versus 



the cluster size n (n = 2 — 10) was found to be in good 
agreement with several previous calculations, in particu- 
lar with ab — initio calculati ons as well as with experi- 
mental data for Cn{n = 2 — 5 ^^ * ^"*^ '. In the present paper 
we extend our previous work and calculated the ground 
state energies for Cn{n = 11 — 55) carbon clusters. A de- 
tailed study will be presented of the possible geometries 
of neutral Cn{n = 11 — 55) carbon clusters and we de- 
termine the ground state, using the energy minimization 
technique with the Brenner potential function using the 
CG method described in our previous work (Ref. 26). 
This method is used to obtain equilibrium geometries, 
bond length, and binding energies for carbon clusters of 
arbitrary size. The number of possible stable isomers in- 
creases very rapidly with increasing cluster size, and it is 
not possible to study all those clusters for n beyond 30 
atoms. For this reason in this work we limit ourselves to 
the lowest energy configurations and pay special atten- 
tion to hexagonal and trigonal shaped graphene fiakes. 
The effect of passivation of the edge atoms by H on the 
energy and stability of the clusters is studied in detail. 

This paper is organized as follows. In Sec. II, a de- 
scription is given of the inter-atomic potential used in 
our calculation. The ground state geometries and ener- 
getics of isomers of small carbon clusters Cn{n = 11 — 55) 
predicted in this study are compared with available the- 
oretical and experimental data of others in Sec. III. The 
effect of H-passivation of the edge atoms on the energy 
and the configuration is investigated in Sec. IV. The sta- 
bility of carbon clusters are discussed in Sec. V. The 
energetics of trigonal and hexagonal shaped clusters are 
studied in Sec. VI. In Sec. VII we summarize our main 
results. 



II. SIMULATION METHOD 

The computational technique and the Brenner second- 
generation reactive empirical bond order (REBO) 
potential^^ function between carbon atoms are similar as 
those used in our previous work (Ref. 26). Here we will 
concentrate on the main differences with the approach 
of Ref. 26. The values for all the parameters used in 
our calculation for the Brenner potential can be found 
in Ref. 28 and are therefore not listed here. The Bren- 
ner potential (REBO) terms in this work were taken as 
follows: 



V^{r) = f{r) Y, B^e-^-^ 



(3) 
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where the cut-off function f^{r) is taken from the switch- 
ing cutoff scheme 
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A^Q^a^Bn^ /3n^{l < n < 3) are parameters for the 
carbon-carbon pair terms. Here n is the type of bonds 
(i.e. single, double and triple bonds). 

The empirical bond order function used in this work is 
written as a sum of terms: 



bij = ^Mr+b'^r]+bi, 



(5) 



where the functions b^~^ and b^~^ depend on the local 
position and bond angles determined from their arrange- 
ment around each atom {i and j, respectively) and is 
governed by the hybridization of the orbitals around the 
atom. Since in the present work we restricted ourselves 
to quasi-planar structures, the conjugate compensation 
term Fij\Ni,Nj,N^p^) was omitted following Ref. 13, 
which is now given by 
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Here the angular function G{cos{Oijk)) modulates the 
contribution of each nearest neighbour and is determined 
by the cosine of the angle of the bonds between the atoms 
i-j-k. 

Some angles between two carbon atoms are less than 
109.47°, as for instance in the case of pentagonal rings 
(Ref. 26). Therefore the revised angular function gc for 
109.47° and 0° (6-7) was used 

gc = Gc{cos{0))^Q{N^)[^c{cos{0))-Gc{cos{0))] (7) 
with 
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Here £'& is the average binding energy in electronic 
volt (eV), V^{rij) and V^{rij) are a repulsive and an 
attractive term, respectively, where r^ is the distance 
between atoms i and j given by 



V'^{r) = f{r){l + Q/r)Ae- 



(n,)] (1) /J ^c < 3 2 

Qi{N^) = < [1 + cos(27r(A^f - 3.2))/2] 3.2 < iVf < 3.7 
[ ATf > 3.7 

(8) 
where N^ is the number of neighbours of carbon atom i. 
^c{cos{0)) is a second order spline which was determined 
for angles less than 109.47° which is coupled to G{cos{0)) 
through the local coordination, and it retains the value of 
(2) first and second derivatives at ^ = 109.47° of the original 



function G{cos{0)). The function b^j in equation (5) is 
well described in Ref. 28. 

We calculated all possible planar structures which are 
close to the ground state, and in particular we will con- 
centrate on graphene-like structures. Moreover we stud- 
ied different graphene structures with various shape and 
topology in order to learn which of them has the low- 
est energy, i.e. is the ground state. Furthermore, H- 
passivation of the edge atoms will also be considered. 



III. CONFIGURATIONS AND ENERGIES OF 
CARBON CLUSTERS 

The optimized geometries are shown in Fig. 1. Table 
I lists the energy of the clusters per atom. Their shapes 
and the corresponding bond lengths between the atoms. 
Some data from the literature are given in parentheses 
(Ref. 13) and in brackets (Ref. 22) for the sake of com- 
parison. 

As can be seen from Fig. 1 the most stable planar clus- 
ters are close-packed structures that have mainly zigzag 
edges. This conclusion is in agreement with previous 
works J ^-*^'^^' We found that for n < 19 the most stable 
structure is a mono-ring which contrasts with previous 
results of Ref. 13 where such a single ring configuration 
was found to be the ground state only up to C12. R.O. 
JoneP^ used the density functional technique and found 
mono-ring structures up to Ciq. M. Sawtarie et. al.^^ in- 
vestigated the relative stabilities of ring, multi-ring and 
cage structures of C20 clusters. Raghavachari et. al.^ 
considered a mono-ring structure (with D2/1 symmetry) 
for C24. C.H. Xu et. al obtained mono-ring structures 
up to CiqW^ D. TomanelP^ obtained, by combining an 
adaptive simulated annealing method and a simple tight- 
binding-type Hamiltonian for the total energy, mono-ring 
structures up to n=22. These different maximum n val- 
ues for which a mono-ring is the ground state is a conse- 
quence of the different theoretical formalisms that have 
been used. It also implies that the accuracy with which 
the sp-hybridization is described is crucial in determining 
the stability of these carbon nanorings for n > 11. 

In Fig. 2(a) we plot the binding energy per atom as 
function of the number of atoms in the cluster and com- 
pare it with other calculations and with available experi- 
mental data. Notice that the experimental and ab-initio 
results are close to our results which exhibit somewhat 
less scatter. For n > 5 our results start to deviate from 
the GA/SA data which are lower in energy. The latter 
is a consequence of the different interatomic potential (i. 
e. earlier generation of Brenner potential) that was used 
(the same remark holds for the results of Ref. 22). We 
notice pronounced oscillations in the GA/SA results for 
n > 15 which are related to the shape of the clusters 
and the position of the polygons inside certain clusters. 
These oscillations are absent in our results and in those 
of Wang et al.^^ This casts doubts on the accuracy of 
the energies found from the genetic algorithm of Ref. 13. 



Fig. 2(b) shows our results (full squares) for the binding 
energy over the full range of investigated cluster sizes. 
We can divide Fig. 2(b) in five zones. In the first zone 
we found single ring clusters Cii_i8 which are structures 
having Dnh {n = 11 — 18) symmetry without variations 
in the bond lengths. The bond length decreases with 
increasing n and lies in the range 1.348-1.338 A, which 
is typical for a double-bond, and which compares with 
the bond lengths 1.334-1.338 A for Cu and 1.333-1.336 
A for C12 obtained by the GA/SA technique. In Ref. 
24 these mono-ring structures with the same symmetry, 
were predicted up to n < 17 with the binding energy of 
Ci7 equal to -6.1668 eV/atom which compares with our 
result -6.0093 eV/atom. 

The reason for the stability of such mono-ring config- 
urations for large n was given by M. Sawtarie et al.^ 
The larger ring diameter of single ring structures results 
in a large angle between adjacent C — C bonds which 
minimizes the overlap of tt electrons in the plane of the 
carbon atoms, that would otherwise have resulted in a 
symmetry lowering through in-plane distortions. 

The second zone starts with a drop in the derivative 
of the binding energy where the ground state structure 
transforms from a single ring to a planar graphene-like 
structure. C19 is a structure with five hexagons, which 
forms zigzag edges with the average binding energy - 
6.2220 eV/atom. There are small variations in the bond 
lengths, i. e. 1.371-1.451 A, which compares with a bond 
length of 1.42 A for bulk graphene. Thus, the transfor- 
mation from a single ring structure to a five hexagon 
graphene structure causes a drop in the binding energy. 

The ground state the C20 cluster consists of a bowl- 
like structure, which agrees with the result of Ref. 13. 
It has one central pentagon and five hexagons surround- 
ing it, with bending |z|^aa:=0.965 A (see Fig. 3). The 
bond lengths have a smaller variation as compared to 
those of Ci9 (Table I), with a maximal difference of 0.061 
A. We found that Cu and C20 carbon clusters can form 
metastable structures with rhombic polygons. In the case 
of C21 the most stable structure is the planar multi-ring 
isomer which is constructed from one central pentagon, 
surrounded by four hexagons and one heptagon. C22 is 
similar to C21 where one hexagon is replaced by a hep- 
tagon. Notice that these heptagons in the structure pre- 
vent the bending of small clusters with a pentagon in 
the centre. Such clusters were obtained by GA/SA for 
C20 (Ref. 13), and by the TQM method for the C21 and 
C22 configurations. The C23 and C24 are small hexagonal 
shaped structures. (^24 consists of purely seven hexagons 
while C23 has a defect, i.e. a pentagon at the edge. The 
extra atom that sits on the top of the armchair edge is 
responsible for an ac(56) edge reconstruction.^^ Thus C23 
has one hexagon less than C24. The bond length differ- 
ences in the structure are 0.093 A for C23 and 0.058 A 
for C24. Since, the reconstructed edge of C23 is symmet- 
ric, the bond lengths between armrest atoms in the edge 
is 1.417 A, and bond length of the C-atom with dan- 
gling bond is 1.469 A. The second zone finishes at the 
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FIG. 1. The ground state configurations for Cn{n = 1 — 55) 
clusters. Tiie underlined numbers are for clusters whose 
configuration changes when H-passivated. Pentagons (hep- 
tagons) are coloured red (blue). 



C24 cluster. C24 is a planar graphene structure and it is 
the evolutionary form of C23 with seven hexagonals and 
a maximum bond length difference of 0.058 A. The C24 
cluster is hexagon shaped with zigzag edges. One can di- 
vide the structure into two types of six C-atom hexagons: 
inner and outer one. The inner polygons bond lengths 
are 1.443 A, the outer ones are in the range 1.399-1.391 
A, the most outer bonds lengths, which involve both sub- 
lattice atoms, are shorter in comparison with the others. 

The bonds that connect the inner ring with the outer 
one have a length of 1.499 A and are the longest bonds 
in this graphene structure. The average binding energy 
of the edge atoms in C24 is -5.8677 eV/atom. 

The binding energy increases at the C25 graphene 
structure. Its ground state consists of seven hexagons as 
C24, but with a different arrangement of polygons. Thus 
C24 is a close packed structure, which purely consists 
of hexagons, while the hexagons in C25 are far removed 
from each other. C25 has two type of edges: four zigzag 
and one armchair edge. It is the first structure which has 
two different type of edges (the Ci9_24 clusters have only 
zigzag edges). The bond length differences in the cluster 
is 0.081 A. 

The planar C26 and C27 clusters are evolutionary forms 
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FIG. 2. The binding energy Eb per atom versus the number 
n of C-atoms in the clusters, (a) Comparison between our re- 
sults and those from GA/SA, experimental data, ab — initio, 
and TQM calculations, (b) The red circles are the bind- 
ing energy of carbon atoms for ground state structures with 
H-passivation, and the open triangles are the corresponding 
(metastable) structures without H-passivation. The ground 
state energy of the pure C-clusters is given by the solid square 
symbols which can be divided in groups as shown by the ver- 
tical dashed lines. The horizontal dashed line separates the 
H-passivated clusters into two groups. 



of C23 and C24 structures. C26 has a defect due to the 
removal of one atom from the zigzag edge of C27, and thus 
it has one pentagon. C27 consists of eight hexagons and it 
is a graphene- type structure. The bond length differences 
for C26 are 0.099 A and 0.08 A for C27. As in the previous 
structures the longest bond length corresponds to the 
outer bonds with average bond length 1.446 A. 

C28 and C29 are the evolutionary forms of the C20 and 




FIG. 3. Dependence of the clusters bending \Z\max (A) 
function of the cluster size n. 



C21 clusters, which are obtained by adding six atoms to 
the edge. We have a cap-like C28 and planar C29 con- 
figurations (see Fig. 3), with the bond length differences 
of 0.075 A and 0.078 A, respectively. The bending, of 
C28 is |z|max=l-535 A. The longest bond lengths appear 
in the pentagon with average length of 1.446 A as in 
C26' The pentagon in C29, like in C21, is surrounded 
by four hexagons and one heptagon. The outer torsion 
bonds of the heptagon have the shortest bond length, i.e. 
1.370 A. The next structures with the number n = 30 - 
33 are planar, evolutionary forms of each others. Their 
basic configuration is C24. All of them are, graphene- 
type clusters, except C31, which exhibits an ac(56) type 
of edge reconstruction. Since the reconstructed edge is 
symmetric as in C23, the bond lengths between armrests 
atoms is 1.418 A and for the dangling bonds it is 1.471 
A. The bond lengths in C30 differ at most by 0.079 A. 
There are two types of edges in the structure: armchair 
and zigzag. The length of the cluster along the arm- 
chair edge is 6.057 A and along the zigzag edge is 5.983 
A. Nevertheless the middle width length of the struc- 
ture is 6.964 A that means, the cluster is bended inside 
by zigzag edges. The bond length differences in C31 is 
0.101 A. The length of the cluster is 7.489 A and the 
widths are 5.265-7.003-5.734 A. This cluster is bended to 
the outside. The average energy per atom for the edge 
with pentagon is -5.7242 eV/atom, and for the edge with 
hexagon is - 5.8751 eV/atom. The bond lengths differ- 
ence in C32 is 0.08 A. The width of the cluster is 5.712 
A and length is 7.397 A, and it is bended toward the 
outside like C31. All edges are of the zigzag type with 
average binding energy -5.92 eV/atom. 
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FIG. 4. Comparison between our results and those from 
GA/SA (Ref. 13) for the number of polygons: (a) pentagon, 
(b) hexagon and (c) heptagon in the cluster. Open triangles 
are those for H- passivated clusters. 



The binding energy exhibits a jump of 0.0311 eV/atom 
at n = 33 were region IV starts. The bond length differ- 
ences in C33 is 0.078 A. C33 has two types of edges: one 
armchair and three zigzag edges. Two of the zigzag edges 
have the same size and they are shorter than the other 
one. The length of the cluster along the armchair edge 
is 5.712 A and along the zigzag edges are 7.417-9.956 A. 
The cluster has a trapezoidal shape. 

To understand better the jump in the binding energy 
per atom when we go from C32 to C33 we show in Fig. 5 
the detailed spatial configuration of the atoms, together 
with the absolute value of the energy of each C-atom in 
the cluster and the distance between the C-atoms. No- 
tice that ternary atoms (i.e. atoms with three carbon 





FIG. 5. Spatial configuration of tfie (a) C32 and (b) C33 
clusters with atomic energies (in eV/atom) and bond length 
(in A) . The same colors are used for atoms and bonds having 
practically the same value for the energy and bond length, 
respectively. 



niegbours) have a larger energy, i.e. E ^ 7.35 eV, which 
depends slightly on the exact positions in the cluster. 
Those atoms are connected with three covalent bonds of 
length 1.44 A. At the edge we can have secondary atoms 
that have a substantial lower energy, typically about 5.13 
eV with shorter bond length, of about 1.40 A. 

The cluster C32 has only zigzag edges with in total 22 
edge atoms of which 14 are secondary atoms (i.e. atoms 
with only two bonds). The other cluster C33 has both 
zigzag and armchair edges with in total 24 edge atoms of 
which 15 are secondary atoms. Thus when we go from 
C32 to C33 the number of edge atoms increases with one. 
It is the latter increase that is responsible for the smaller 
average binding energy per atom of the cluster. 

Here, a small comment is in order about the two 
types of evolutionary forms of small clusters that we 
have found. For instance, the ground state configuration 
of the clusters Cn-is, C20-22, ^23, 24, C'26-33, ^^35-37, 
C'40-42, £^46, 47, and C49,5o are obtained by adding only 
one atom to a certain ground state configuration. The 
other type of configurations result from an increase of 
polygons, such as the structures C20, ^^23, C26, C'34, C38, 
and C44. Thus, we can distinguish them as atomic and 
polygonal evolutions. One of those evolutions leads to 
the creation of planar structures, when the pentagon as 
a defect is inserted at the edge resulting in an edge re- 
construction. The other type of evolution leads to the 
creation of bended structures, when the pentagon is in 
the centre of the clusters (see Figs. 1 and 3). So C34, 
in addition to C26 has additional three hexagons. The 
structures bond length differences of 0.11 A. This large 
difference is due to the pentagon in the cluster with bond 
length 1.444 A. All edges in C34 are of zigzag type. 

The next clusters, i.e. (735-37, are atomic evolutionary 
forms and are all planar configurations. Two of them, 
C35 and C37 consist of 11 and 12 hexagons, respectively, 
with bond length differences of 0.081 A for C35 and 0.08 
A for C37 as in C32. C35 has two types of edges: arm- 
chair and zigzag. C37 has only zigzag edges. The average 
binding energy £^5 of the edge atoms in the cluster is - 
5.9380 eV/atom. Csq in addition to its eleven hexagons 
has one pentagon with an armchair ac(56) type of edge 
reconstruction. The difference of bond lengths in this 
structure is 0.103 A. This large difference between bonds 
is due to the pentagon in the cluster (i.e. as in C35) 
which has an average bond length of 1.453 A. The Css 
structure is the polygonal evolutionary form of C20, with 
additional seven hexagons. The bond length differences 
in this structure is of 0.078 A. All edges are of the zigzag 
type, and the structure is bended by |z|maa;=l-891 A. C39 
structure is a unique cluster with one heptagon, having a 
planar structure with bond length differences of 0.087 A 
and only zigzag edges. The C40 structure is a planar one, 
and it is a polygonal evolutionary form of C37 with one 
armchair edge. Thus, C40 has two types of edges. The 
bond length differences is 0.081 A. The two planar con- 
figurations C40 and C41 are atomic evolutionary forms 
of each other. C41 has one pentagon on the edge, and 



thus has only zigzag edges. The bond length differences 
are 0.104 A for C40, and 0.079 A for C41. Adding one 
atom to the pentagon of the reconstructed ac(56) edge 
of C41 we obtain a close packed structure as the ground 
state having only zigzag edges (see Fig.l). The struc- 
tures bond length difference is 0.079 A. The fourth zone 
finishes with the C42 cluster, where the average binding 
energy of the edge atoms is -5.9538 eV/atom. 

The C43 structure is a planar graphene-type cluster 
with bond lengths difference of 0.079 A as in the C42 
structure. Note that the jump in the binding energy per 
atom is 0.168 eV/atom which is due to the fact that C43 
has two types of edges because the average binding en- 
ergy per atom in the armchair edge is larger than in the 
zigzag edge. The reason is that armchair termination in- 
volves both types of atoms (i. e. both sub-lattice sites). 
Thus, two types of atoms are close to each other and form 
pairs. The average bond length in this case is 1.373 A, 
which is shorter than the average bond length for zigzag 
edges (1.41 A). The C44 structure is a polygonal evo- 
lutionary form of C20' The bond lengths have a differ- 
ence of 0.076 A and the cluster is strongly bended with 
|^|max= 3.290 A (see Fig. 3). All edges are of zigzag type. 
C45 is a planar graphene-type cluster which is a polyg- 
onal evolutionary form of C40 with the addition of two 
hexagons. The bond lengths have a difference of 0.079 A, 
as in the C42 structure, and has two types of edges. The 
planar C46, C47, C49, and C50 configurations are atomic 
evolutionary forms. C46 and C49 have an edge defect as 
pentagon, but with different reconstructions: the edge 
defect of C46 is of ac(56) type while the pentagonal de- 
fect in C49 was formed by removing one atom from the 
hexagon. The bond length differences are 0.106 A, and 
0.081 A respectively for C4Q and C47, and 0.079 A for 
both C49 and C50. 

The ground state configurations for n > 50 are 
graphene-like clusters. However, C51 and C53 have a de- 
fect on the edge. The defect of the C53 cluster is formed 
due to the lack of one C atom in the hexagon which was 
positioned in the centre of the zigzag edge. In all cases 
such kind of ac(56) or ac(67) defects do not lead to a 
bending of the cluster. The C52, C54, and C55 structures 
are purely graphene-type clusters. C52 and C54 consists 
of zigzag edges, and C54 is a hexagon shaped configura- 
tion. C55 looks like a tiny nanoribbon with armchair and 
zigzag edges. 

Each time the cluster possesses both types of edges 
the binding energy per atom jumps up (see Fig. 3 (b)). 
Armchair edges change the cluster structure, and make 
them extended, rather then close-packed. As shown in 
Fig. 2(b) and Table I the binding energy exhibits a drop 
when the clusters are close-packed, having only zigzag 
edges. Moreover armchair edges cause an increase of the 
number of secondary carbon atoms (see Fig. 8). Thus, 
when we have a larger number of atoms with only two 
neighbours it results in a smaller binding energy. We 
found that 42% of the lowest energy configurations are 
graphene-type structures. 



It is clear from Fig. 4, that the number or the type of 
polygons in the clusters can not be the cause of oscilla- 
tions in the binding energy. The oscillations are due to 
the position of the polygons in the structure, and whether 
the structure has zigzag or armchair edges. 

9 % of the investigated clusters are bended structures 
and in all of them there are pentagons in the centre (see 
Figs. 1 and 3). Those pentagons are surrounded by 
hexagons and bend them such that their bond length cor- 
responds to a double bond. The smallest bended cluster 
in the range n= 11-55 is the C20 cluster and the largest 
one is C44 (see Fig. 3). All the clusters after C44 are pla- 
nar configurations. Notice that the buckled structures do 
not affect the binding energy oscillation. 

In Figs. 4 (a, b, c) we plot the number of differ- 
ent polygons, i.e. pentagons (P(5)), hexagons (P(6)), 
and heptagons (P(7)), constituting the cluster as func- 
tion of the size of the cluster. Notice that the number 
of hexagons steadily increases while the number of pen- 
tagons and heptagons are mostly one and occasionally 
two. In the theoretical results of Ref. 13 the pentagons 
and heptagons play a more important role as compared 
to our work. The number of pentagons in our results are 
substantially smaller than in the results of Ref. 13, ex- 
cept for the C44 cluster in our work, and except C21 in the 
GA/SA work (both are buckled). However, in our results, 
starting from C22 it is clear that the hexagons dominate 
which at the end will result in graphene-like clusters (see 
Fig. 4(b)). We found that clusters which mainly consist 
of pentagons and heptagons are energetically not favor- 
able, i. e. they are metastable configurations. The sin- 
gle heptagon ring is most stable with average energy per 
atom -5.5081 eV/atom, which compares with the single 
pentagon with average binding energy -4.742 eV/atom 
(see Ref. 26). However in non mono-ring structures the 
heptagons are positioned at the outer shells while the 
pentagons are located either in the center of the cluster, 
or at the outer shell. The reason is that the heptagon 
has a smaller bond length (1.352 A) which is close to 
the double bond length, which compares with pentagons 
which have 1.437 A (see Ref. 26). Both, our results and 
those from the GA/SA work find that hexagon rings play 
a dominant role in the ground state of clusters. 



IV. EFFECT OF H PASSIVATION OF THE 
EDGE ATOMS 

Often edge atoms are passivated by hydrogen atoms, 
which will remove the dangling bonds. To our knowledge, 
the effect of H atoms on the structure of carbon clusters 
has not been studied yet. In what follows we perform a 
systematic study of the influence of H passivation of the 
edge atoms on the stability of carbon clusters. 

We attach H atoms on the edge C atoms of all ground 
state, and metastable configurations which we studied in 
previous section. From Fig. 2(b) and Table I we see that 
the binding energy £^5 of hydrogen passivated clusters is 



TABLE I. The binding energy per atom as function of the number of atoms in the (C-C) and (C-H) clusters (the numbers 
between parentheses are from Ref. 13, and those between brackets are from Ref. 24). The cluster shape is labeled SR for 
single rings, G for graphene, Gcp for close packed, DR for double ring, MR for multi-ring, B for bended structures. The fourth 
column gives the range of bond lengths appearing in the cluster (results between parentheses are those from Ref. 13). The last 
three columns are the equivalent results for the H-passivated carbon clusters. 



n 


E^- 




Shape 


Dcc 


E^- 


Shape 


Dcn 




(eV/atom) 


(C-C) 


(A) 


(eV/atom) 


(C-H) 


(A) 


11 


-5.8647 




SR 


1.348 


-7.5911 


DR 


1.377-1.428 


12 


-5.9041 




SR 


1.346 


-7.5857 


DR 


1.384-1.415 


13 


-5.9341 




SR 


1.344 


-7.5557 


Gcp 


1.382-1.428 


14 


-5.9595 




SR 


1.342 


-7.5606 


G 


1.383-1.430 


15 


-5.9794 




SR 


1.341 


-7.5671 


MR 


1.373-1.444 


16 


-5.9957 




SR 


1.340 


-7.5406 


Gcp 


1.381-1.427 


17 


-6.0093 




SR 


1.339 


-7.5382 


MR 


1.372-1.450 


18 


-6.0207 




SR 


1.338 


-7.5498 


G 


1.383-1.430 


19 


-6.0681 




Gcp 


1.371-1.451 


-7.5299 


G 


1.381-1.427 


20 


-6.1041 ( 


;-6.475) 


MR 


1.390-1.451 


-7.5282 


MR 


1.371-1.449 


21 


-6.1424 1 


-6.44] 


MR 


1.367-1.457 


-7.4361 


MR 


1.383-1.443 


22 


-6.1674 1 


-6.465] 


MR 


1.360-1.458 


-7.5429 


G 


1.380-1.427 


23 


-6.1992 1 


-6.49] 


MR 


1.376-1.469 


-7.5265 


G 


1.381-1.429 


24 


-6.2368 1 


-6.51] 


Gcp 


1.391-1.449 


-7.5121 


Gcp 


1.381-1.426 


25 


-6.2013 




B 


1.370-1.451 


-7.5162 


G 


1.381-1.428 


26 


-6.2533 1 


-6.525] 


MR 


1.378-1.477 


-7.5200 


G 


1.380-1.428 


27 


-6.2727 1 


-6.55] 


G 


1.369-1.449 


-7.5075 


G 


1.381-1.426 


28 


-6.2834 1 


-6.58] 


B 


1.376-1.451 


-7.5116 


G 


1.380-1.427 


29 


-6.3083 1 


-6.65] 


MR 


1.370-1.448 


-7.4491 


MR 


1.363-1.489 


30 


-6.3148 




G 


1.369-1.448 


-7.5039 


G 


1.381-1.426 


31 


-6.3429 




MR 


1.370-1.471 


-7.4390 


B 


1.368-1.457 


32 


-6.3684 




Gcp 


1.368-1.448 


-7.4974 


Gcp 


1.381-1.426 


33 


-6.3373 




G 


1.369-1.447 


-7.5010 


G 


1.381-1.428 


34 


-6.3704 




MR 


1.372-1.482 


-7.4302 


MR 


1.383-1.462 


35 


-6.3963 




G 


1.367-1.448 


-7.4952 


G 


1.381-1.426 


36 


-6.4180 




MR 


1.371-1.474 


-7.4280 


MR 


1.383-1.463 


37 


-6.4380 




Gcp 


1.368-1.448 


-7.4899 


Gcp 


1.381-1.426 


38 


-6.4365 




B 


1.374-1.452 


-7.4933 


G 


1.381-1.426 


39 


-6.4473 




MR 


1.370-1.457 


-7.4311 


MR 


1.383-1.470 


40 


-6.4575 




G 


1.367-1.448 


-7.4885 


G 


1.381-1.426 


41 


-6.4748 




MR 


1.370-1.474 


-7.4294 


MR 


1.383-1.464 


42 


-6.4913 




Gcp 


1.368-1.447 


-7.4943 


G 


1.381-1.476 


43 


-6.4745 




G 


1.368-1.447 


-7.4875 


G 


1.381-1.426 


44 


-6.4860 




B 


1.378-1.454 


-7.4905 


G 


1.381-1.426 


45 


-6.5054 




G 


1.368-1.447 


-7.4833 


G 


1.381-1.426 


46 


-6.5195 




MR 


1.368-1.474 


-7.4864 


G 


1.381-1.426 


47 


-6.5336 




Gcp 


1.367-1.448 


-7.4797 


Gcp 


1.381-1.426 


48 


-6.5263 




G 


1.368-1.447 


-7.4828 


G 


1.381-1.426 


49 


-6.5393 




MR 


1.368-1.447 


-7.4420 


MR 


1.381-1.426 


50 


-6.5520 




G 


1.368-1.447 


-7.4792 


G 


1.381-1.426 


51 


-6.5650 




MR 


1.369-1.474 


-7.4326 


MR 


1.383-1.463 


52 


-6.5676 




Gcp 


1.368-1.447 


-7.4761 


G 


1.381-1.426 


53 


-6.5918 




MR 


1.372-1.497 


-7.4289 


MR 


1.384-1.469 


54 


-6.5973 




Gcp 


1.367-1.439 


-7.4815 


G 


1.381-1.426 


55 


-6.5753 




G 


1.368-1.447 


-7.4758 


G 


1.381-1.427 



larger and decreases with increasing n, while the bind- 
ing energy of the pure C-clusters is a decreasing function 
of n. Both curves converge to Eq = —7.35 eV/atom in 
the limit n ^ oo. We found that our results for the 
pure C clusters could be fitted to E = Eo-\- ac/n'^'^ with 
Qfc=5.68 eV/atom and for the H-passivated clusters with 
ac= - 0.77 eV/atom, where we did the fitting over the 
range 17 < n < 56. The square root dependence is a 
consequence of the fact that the number of edge atoms 
is proportional to n^-^. The binding energy exhibits os- 
cillations as function of n around these fitted curves. 

The horizontal dashed line in Fig. 2(b) divides the H 
passivated clusters into two regions. Below this line we 
have graphene type clusters (see Fig. 6) with the excep- 
tion of Cii, Ci5, Ci7 and C20 which all have one hep- 
tagon. Above this line those clusters have a defect at the 
edge, i. e. one pentagon. Notice that the ground state 
configuration of the small clusters consists no longer of 
a mono-ring but are mostly graphene-like clusters. We 
found double, triple, and quarted clusters of hexagons. In 
the case of pure carbon cluster calculations these struc- 
tures were metastable configurations (see Fig. 2). H pas- 
sivation turned them into ground state configurations. 

The ground state configurations with hydrogen pas- 
sivated edge atoms are depicted in Fig. 6. The shape 
together with the shortest and the longest bond length 
as given in Table I. Notice that the H passivated clus- 
ters have a smaller range of bond lengths as compared to 
the clusters without H passivation. Similar as on previ- 
ous section, we have two kinds of evolutionary forms of 
structures with increasing size. 

The non-graphene-like structures can be dividing into 
two types by the position of the pentagon defect. In the 
first type the defect is formed by the lack of a C atom 
at the edge hexagon. That kind of defect appears in all 
small C clusters, resulting in a lower binding energy (see 
Fig. 2). The other type of clusters, have a pentagon 
defect formed by the diffusion of carbon atoms from dis- 
tant armrests to seat positions, thus leading to an ac(56) 
reconstruction of an armchair edge (Ref. 40). 

The largest cluster we studied, i. e. n=55, has the 
same configuration with and without H passivation and 
has a graphene structure. The difference between the 
longest and the shortest bond lengths is 0.079 A for the 
pure C-cluster and 0.046 A for the H passivated cluster. 
The clusters with H passivation are slightly compressed 
as compared to the pure carbon clusters. 

69% of the hydrogen passivated structures are 
graphene- type clusters, and only one ground state con- 
figuration C31 is a bended cluster. 



V. MAGIC CARBON CLUSTERS 

In order to find the magic clusters we calculate the 
second finite difference of the binding energy per atom, 
A2E{n) = E{n - 1) + E{n + 1) - 2E{n), which are 
shown in Figs. 7(a,b). The positive peaks indicate par- 
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FIG. 6. The ground state configurations for H passivated 
Cnin = 11 — 55) clusters. The underlined numbers are clusters 
whose configuration is different from the one of the pure C 
clusters that are shown in Fig. 1. Pentagons (heptagons) are 
colored red (blue). 



ticularly stable clusters compared to their neighbouring 
structures. For both C-C and C-H clusters there is no re- 
lation between the stability of the cluster and their even- 
or odd numbered constructions, as was found in the case 
of 3D structures (Ref. 14). In Fig. 8(a) we can see that 
the small size single ring structures are most stable, and 
the biggest single ring C18 structure is the least stable 
one. An analogous trend for single ring structures was 
found in Ref. 24. The graphene-type hexagonal shaped 
C24 cluster is clearly a magic cluster. Note that C24, C32, 
C42, and C54 configurations are polygonal evolutionary 
forms of the two hexagonal shaped C24 and C54 clusters. 
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FIG. 7. The second finite difference of the lowest energy 
A2E{n) per atom vs. cluster size n. Peaks correspond to 
the most stable structures compared to adjacent sizes, (a) is 
for pure C-clusters and (b) is for C-clusters with H-passivated 
edges. 



FIG. 8. The number of zigzag and armchair edges in the 
ground state as function of the size n of the cluster, (a) is for 
pure C-clusters and (b) is for C-clusters with H-passivated 
edges. 



The highly symmetric clusters are the most stable one. 

Although C25, C33, C43, and C48, are graphene-type 
clusters they have armchair edges, and are the least sta- 
ble structures. Despite the fact that C52 has only zigzag 
edges, it has one extra edge atoms in contrast to its neigh- 
bours. Thus, this additional atom increases the number 
of secondary carbon atoms and decreases its stability. 

In the case of H-terminated clusters we see from Fig. 
7(b) that C22, C'ao, C35, and C52 are very stable clusters, 
i. e. magic number clusters, and are polygonal evolution- 
ary forms of the C22 and C52 configurations. We can see 
that the least stable clusters are those having a pentagon. 

In order to study the effect of the type of edges on 
the stability of the clusters we compare the number of 
zigzag edges in a ground state structure with the number 
of armchair edges. We found that zigzag edges occur 
more frequently in both type of ground state clusters 
(see Figs. 8(a-b)). Some structures exhibit solely zigzag 
edges, while no ground state clusters where found with 
only armchair edges (see Fig. 8). 

We can see that each cluster has on the average 5 zigzag 



edges and for n < 43 they have 1 or armchair edges for 
the pure C clusters and a maximum of 3 armchair edges 
for the H passivated clusters. The number of secondary 
carbon atoms for both type of clusters versus cluster size 
n are shown in Fig. 9. The number of secondary carbon 
atoms in the pure C clusters in C18, C25, C33, and C43 
are more than their neighbours. This leads to an increase 
of the number of atoms with higher energies. 

Thus, the abundance of secondary carbon atoms is the 
cause of the structure weakness. 

Generally, H passivated clusters have more secondary 
carbon atoms resulting in an increase of the number of 
ternary carbon atoms. Those atoms have a lower bind- 
ing energy as compared to the case of pure C clusters. In 
all cases when the clusters have a smaller number of sec- 
ondary carbon atoms, they have a larger binding energy 
(see Fig. 2). 
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FIG. 10. Nanodisks with zigzag edges for trigonal and with 
zigzag and/or armchair edges for hexagonal shaped struc- 
tures. 



FIG. 9. The number of secondary carbon atoms in the ground 
state as function of the cluster size n. 



VI. TRIGONAL AND HEXAGONAL SHAPED 
NANODISKS 

Highly symmetric nanodisks have been recently the 
focus of theoretical research^^ ^^. Such nanodisks are 
nanometer scale disk-like graphene materials which have 
closed edges. In particular the electronic structure of 
nanodisks that have the same type of edges, i.e. zigzag 
or armchair, have been investigated. Here we limit our- 
selves to trigonal shaped nanodisks that have zigzag 
edges^ and hexagonal shaped structures^ that can have 
either zigzag or armchair edges (see Fig. 10). The 
interest in those nanodisks originates from the predic- 
tion that zigzag terminated nanodisks have degenerate 
zero-energy electronic stateP^l^ that leads to metal- 
lic ferromagnetism^^ ^^. Recently, STM measurements 
have given evidence for the existence of metallic zero-gap 
graphene flakes with zigzag termination^^. 

We show here that these interesting and promising 
graphene nano-islands are not necessarily the ground 
state configuration for the given number of C-atoms and 
in such a case will be metastable. Nevertheless, it is 
possible to realize such stable and metastable graphene 
nanostructures with zigzag or armchair edges through 
Joule heating of nanoribbons^^, electron bombardment^^ 
or by crystallographic specific etching techniques^ of a 
graphene layer. 

In the theoretical investigation of nanodisks one as- 
sumes that they are build up with hexagonal plaquettes. 
This does not necessarily lead to the most energetic fa- 
vorable state and sometimes the insertion of pentagons 
and/or heptagons may lower their energy. A typical ex- 
ample is the buckyball. We showed already that it is not 



necessary to have a three-dimensional structure for the 
occurrence of non-hexagonal polygons and that even a 
perfect planar structure can have such polygons. 

We consider the trigonal and the two hexagonal shaped 
clusters depicted in Fig. 10. The trigonal structure has 
zigzag edges with a total number of n = l-\-4:Nside-\- N^id^ 
carbon atoms where Ngide is the number of hexagonal 
plaquettes along each side of the triangle. For Nside='^i 
2, ..., 5 we have n=6^ 13, 22, 33, 46 respectively. The 
first one, i.e. n=6, is a single hexagon which was found 
to be the ground state (see Table I). In Table II we also 
show the energy of the unrelaxed cluster where the C-C 
distances are taken equal to the bulk value, a=1.42 A. 
Notice that the relaxation of the configuration, in par- 
ticular of the edge atoms, lowers the energy with about 
0.4%. All trigonal shaped structures shown in Fig. 10 
with the exception of n=6 are metastable state configu- 
rations (see Table II). 

The hexagonal nano-island with zigzag edges has n — 
QN^ize carbon atoms. For Nside='^^ 2, 3 we have respec- 
tively n=6^ 24, 54. We found that all these 3 configura- 
tions are ground state structures. This type of clusters 
are perfectly close-packed planar structures with purely 
zigzag edges. 

The hexagonal nano-island with armchair edges has 
n = 6{l — 3Nside-\-^J^side) carbou atoms. The C42 struc- 
ture is a metastable state cluster, while the ground state 
is a close-packed graphene- type structure (see Fig. 1). 

As we have noted above, the trigonal structures have 
zigzag edges, and they are metastable configurations 
when they consist of only carbon atoms. We studied 
those clusters with H atoms attached to their edges, 
and found that Ce, and the C13 structure are ground 
state, and (722, Css and C46 H passivated structures are 
metastable clusters. 
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TABLE II. The binding energy per atom (in eV/atom) as 
function of the number of sides in the cluster (N) and n is 
the total number of C-atoms. We compare the ground state 
(GS) energy per atom of the pure carbon clusters with the 
H-passivated clusters. For the pure carbon clusters we also 
present the results for the unrelaxed situation (nr). 



Nsize n 


EV 


Eb 


Ei,{GS) 


Eb 


Eb{GS) 




(C-C) 


(C-C) 


(C-C) 


(C-H) 


(C-H) 


Trigonal 












1 6 


-5.2682 


-5.2869 


-5.2869 


-7.6233 


-7.6233 


2 13 


-5.8115 


-5.8427 


-5.9349 


-7.5557 


-7.5557 


3 22 


-6.1036 


-6.1339 


-6.1674 


-7.5219 


-7.5429 


4 33 


-6.2969 


-6.3251 


-6.3373 


-7.5009 


-7.5010 


5 46 


-6.4371 


-6.4628 


-6.5195 


-7.4862 


-7.4864 


Hexagonal (zigzag 


edges) 








1 6 


-5.2682 


-5.2869 


-5.2869 


-7.6233 


-7.6233 


2 24 


-6.0382 


-6.2368 


-6.2368 


-7.5121 


-7.5121 


3 54 


-6.4534 


-6.5973 


-6.5973 


-7.4730 


-7.4815 


Hexagonal (armchair edges) 








1 6 


-5.2682 


-5.2869 


-5.2869 


-7.6233 


-7.6233 


2 42 


-6.2737 


-6.4628 


-6.4913 


-7.4943 


-7.4943 



The trigonal shaped structures have Ns = (6 + 
ANside) — Nside secondary carbon atoms. The (C — (7)22, 
(C — H)22^ iC — C)33 and {C — H)s3 clusters have the 
same number of secondary carbon atoms A^5=12 and 
A^s=15 respectively. {C — C)i3 is a single ring struc- 
ture, while {C — H)is is a graphene-type cluster. In the 
{C — H)is structure all Ns=9 secondary carbon atoms 
after H passivation of the edge atoms become ternary 
atoms. W e observed, that the small {C — H)2a hexago- 
nal shaped clusters with zigzag edges is the ground state 
structure, while (C — i^)54 becomes a metastable config- 
uration. (C — H)^2 is the ground state structure because 
of the abundance of secondary carbon atoms. 



VII. CONCLUSIONS 

Atomic simulations where used with the Brenner sec- 
ond generation reactive bond order (REBO) inter-atomic 
potential function to study small neutral carbon clusters 
Cn{^^ < ^ < 55). In the present work we are interested 
in quasi-planar structures and therefore omitted the con- 
jugate compensation term Fij{Ni^ Nj^ N^^^^) of the orig- 
inal potential function. We concentrated on the ground 
state configurations. Our results show that Cn is a single 
ring for n up to 18, and multi-rings for 19 < n < 55. We 
found that C20, C28, C38, and C44 are buckled clusters 



with a pentagon in the center. 

The binding energy of the clusters does not exhibit 
even-odd number oscillations. But we found smaller en- 
ergy oscillations which are due to the occurrence of spe- 
cial type of structures, and in particular the kind of edge 
they have (i.e. zigzag or armchair type of edges). It also 
depends on the position of some of the polygons in the 
cluster. If the polygons are far removed from each other, 
then the cluster has a higher energy, and if the polygons 
are closer to each other the cluster has a lower binding 
energy. Similarly, we found that clusters with mostly 
zigzag edges have lower binding energy, especially if they 
are close- packed structures than those consisting of both: 
zigzag and armchair edges. 

The hexagons play an essential role in the formation of 
graphene-like clusters. Pentagons can be in the center of 
the cluster, as well as at the edge, while heptagons were 
never found at the center of the ground state structure. 
Clusters having the graphene structure are not bended 
or rippled but are perfectly fiat. 

Our results (see Fig. 2(a)) for e. g. the energy, are dif- 
ferent from previous works which we can trace back to the 
use of a different inter-atomic potential function. Wang 
et al.l^ and Zhang et al.l^ used the previous version of 
the Brenner inter- atomic potential function which differs 
from ours in the fact that our (REBO) potential function 
contains improved analytic functions and was based on 
an extended database as compared to the early version 
(Refs. 28, 42). The bond energies, lengths, and force con- 
stants for hydrocarbon molecules are significantly better 
described in the second version of the potential. 

We found that trigonal shaped clusters having zigzag 
edges are metastable structures. The hexagonal shaped 
clusters that have only zigzag type edges are close- 
packed ground state structures and are always the ground 
state (except for H passivated clusters), while the same 
shaped clusters that have only armchair type edges are 
metastable structures (see Table II). These conclusions 
are valid except for the trivial case n = 6. 

Our calculations show that some of the metastable car- 
bon clusters can become the ground state configuration 
if H-atoms are attached to the edge atoms. The num- 
ber of pentagons influence the binding energy of those H 
passivated clusters. Clusters with pentagons have higher 
energy. The binding energy of H passivated clusters is 
larger in comparison to the pure carbon clusters but both 
converge to the graphene result in the limit n ^ 00. 

For n > 11 graphene-type of clusters can be formed as 
ground state. This is even more so for n > 41. 
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